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Abstract — The paper studies the dynamics of a conventional 
positive going zero crossing type digital phase locked loop 
(ZC1-DPLL) taking non-ideal responses of the loop constituent 
blocks into account. The finite width of the sampling pulses 
and the finite propagation delay of the loop subsystems are 
properly modeled mathematically and the system dynamics is 
found to change because of their influence considered 
separately. However, when these two are taken simultaneously, 
the system dynamics can be made nearly equivalent to that of 
the ideal system. Through an extensive numerical simulation 
a set of optimum parameters to overcome design limitations 
have been obtained. 

Index Terms— ZC1-DPLL, Delay, DCO Pulse Width, Bifurcation 
Diagram, Layapunov Exponent, Stability Zone, Convergence 
Time 

I. Introduction 

Among different structures of digital phase locked loops 
(DPLLs) positive going zero crossing sampling type loops 
(ZCl-DPLLs) are widely used in practical applications [1,2,3,4]. 
In this DPLL, the input sinusoidal signal is sampled at the 
transition instants of the loop digitally controlled oscillator 
(DCO) signal and the DCO frequency is controlled by the 
filtered version of the sampled input signal. In the 
mathematical model of the ZC1-DPLL one makes a few 
approximations that are valid for an ideal system. For example, 
the width of the sampling pulses is taken to be infinitesimally 
narrow and as such the sampled voltages are considered as 
the instantaneous values at the sampling instants. Further, 
since the next sampling instant of the input signal is to be 
determined by the value of the sampled signal at the present 
instant, some processing has to be done on the sampled 
signal using physical circuits (viz, analog to digital converters, 
digital filters, DCO, etc.). Hence there would be a finite time 
delay in system response [5]. In mathematical modeling of 
the DPLL, this time delay is not considered. 

In this paper our aim is to study the dynamics of a practical 
ZC1-DPLL without making these two approximations. Thus 
we treat the loop as a 'non-ideal' one. The finite width of the 
sampling pulse has been modeled by taking time averaged 
values of sampled signals in the duration of the sampled 
pulses and the finite propagation delay of the closed loop 
circuit is modeled by the inclusion of a pure phase delay 
network in the functional model of the loop. After deriving 
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the system equation incorporating the abovementioned 
considerations, the dynamics of the practical ZC1-DPLL has 
been studied both numerically and quasi-analytically. Suitable 
design parameters have been assigned to take into account 

of the finite width of the sampling pulses (At) an d nonzero 
loop time delay ( r ). It has been observed that when the 
magnitudes of these become more than some critical values, 
the system dynamics becomes chaotic. Further these two 
effects have mutually opposite influence on the loop 
dynamics. From a different angle of view, we can point out 
that a controlled application of loop time delay or non-ideal 
finite width sampling pulses can provide us with predictable 
chaotic signal generators. These systems are of importance 
in chaos-based communication system design. 

II. Derivation of system equation for non-ideal ZC1-DPLL 

The functional block diagram of an ideal ZC1-DPLL is 
shown in Fig. 1 . It comprises of a sampler, a loop digital filter 
(LDF) and a digitally controlled oscillator (DCO). 
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Fig. 1 Functional block diagram of ZC1-DPLL 

For a non-ideal ZC 1 -DPLL we consider the width of sampling 
pulses as At and an over all loop propagation delay r . The 
effect of the finite width of the sampling pulse would be felt 
at the output of the sampler. If the input signal to the sampler 
be, 

x(t) = Aq sin(<y ; f + Q ) ( 1 ) 

Where the noise free input signal x(t) to have an angular 
frequency o, , amplitude Aq and constant phase 6 . The 
sampled signal at the k th sampling instant t = t(k) 
(k = 0,1,2 ) is written as x(k) for an ideal DPLL as, 
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x(k) = x(t(k)) = A) sin( £»,■?(£) + O ) (2) 
However, for finite width sampling pulses the sampler output 
at the k-th sampled instant (SI) would be taken as a time 

average of x(t) around t(k) for a period At ■ Thus, sampled 
signal for a non-ideal DPLL would be given as x(k) where 

f(*)+Af 



Input signal 



x(k) = — f Aq sin( af)dt 

At Hk) 



(3) 



Normalizing the input signal frequency CO i in terms of the 
DCO nominal frequency w as, to i = %co , and writing At in 
terms of the DCO nominal period T , At = rT = r(2x/<D ) . One 



can write, 



x(k) = GiAq s'm[a>it(k) + fytr\ 



(4) 



_ sin(^OT-) 

Here we put, - — — . It is to be noted that the quantity 

Gj accounts for the finite width of the sampling pulses 

through the parameter r(= At/T Q ). As y tends to zero would 
tend to one making the DPLL an ideal one. 

Further it has been already mentioned that to analytically 
model an ideal ZC1_DPLL, the responses of the subsystems 
of a ZC1 -DPLL (viz. sampler and quantizer, LDF and a DCO) 
have been considered to be instantaneous. But in practice, 
they will introduce a time delay, whatever small, in transmitting 
the signals through them and as such in the close loop 
response of the system one has to take into account an 
additional inherent time delay. This delay would definitely 
influence the dynamics of the loop. To derive the system 
equation of a delayed ZC1-DPLL the overall time delay has 
to be modeled by a suitable Z-domain transfer function. In 
the literature, a pure time delay (x) is modeled in analog s- 
domain by the operator exp(-sx), (s = jco). For a small delay, it 
can be approximated as, 



2 



1 + ^ 
2 



(5) 



Performing a bilinear transformation one can get the Z domain 
equivalent of the pure delay given in [5] as, 



1+ az 



(6) 



Where a is a quantitative measure of the time delay, normalized 
to the sampling period of the digital system. Theoretically 
the parameter a can take values in the range o < a < 1 , but for 
simplified derivation of the system equation we may take the 
range as < a < 1 , where a= indicates a delay of one 
sampling period and a = 1 indicates no delay. As expected, 
the delay network thus modeled gives unit amplitude gain 
and only pure phase shift, depending on the magnitude of a. 
The functional model of a "non-ideal" DPLL considered in 
this paper is shown in Fig.2. Here an average block (operating 
for an interval of sampling pulse width At ) and a pure delay 
network have been included. 
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DOr.01.IJCOM.3.1.2 



Sampler 


— » 


Averager 





DCO 




Delay 






4 — 


4 



Fig. 2 Functional block diagram of delayed ZC1-DPLL 

Now, the unmodulated and noise free input signal x(t) to 

the loop is considered to have an angular frequency C0 i , 

amplitude Aq and constant phase 9 Q . The signal is written 

in terms of the DCO nominal frequency C0 Q as, 

x(t) = A Q sm(co Q t + 0(f)) , where 9(t) is defined as 

(co j — O) )t + . x(t) is sampled at the instants determined 
by the DCO signal's positive going transition from one level 
to the other. The output of the LDF at the SI denoted by 

y(k) is used to control the period of the DCO according to 
the algorithm, 

T(k+\) = T(k)-y(k) 

Here T (= 2n I co ) is the nominal period of the DCO and 

T(k+l) is defined as the time elapsed between fc th and 

(k + l)' h SI, i.e. T(k + 1) = t(k + 1) - t(k) . Considering t (0) = , 
one can write, 

t(k) = kT-Y,y(i) (7) 

o 

Now, considering the LDF as a gain element of magnitude 
G and taking the pure time delay network into account one 
can write: 



y(k)- 



a + z 



l + az 



G Q x(k) 



(8) 



Here x(k) stands for sampler output for an ideal ZC 1 -DPLL. 
The phase error between the input signal phase and the loop 

DCO phase is given by (p(k) which is the signal angle at the 

sampling instant t(k) . Thus, 

<p(k) = (a i t(k) + e 

Using (7), one gets (p(k) 

<p(k) = afk) + O =% 0Jat(k) + 6 



--^ka) T -^a Y J y(i) + % 

<=o 



(9) 



49 



Where, £ = (»/»„ ) . For an ideal ZC 1 -DPLL, the sampler out- 
put x(k) = Aq sin (p(k) , Now, after simple derivation of eqn.(9) 
using eqn.(8) one gets the phase governing equation as 
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cp(k + 2) = (1 -aMk + 1) + aqtji) + 2x(l+ a)(£ - 1) - 
£ K (a sin <p(k + 1) + sin q>(k)) (10) 
Here f and K have been substituted in place of 
(ft), /a) ) and AqCOqGq . 

Again due to consideration of finite width of sampling pulses 
the sampler output x(fc)must be replaced by x(k)in system 

equation derivation and after a simple algebra the system 
equation of the non-ideal system may be written as 

<p(k + 2) = (1 - a,ri)<p(k + 1) + anqik) + 2x(l + an) - 1) - 
^K [asii(p{k + 1) + 7i$r) + {1 + (n- l)a 2 Jsii^A:) + ^r)) (11) 
For a non -ideal DPLL, n = 1 and for ideal DPLL a=l, r = 
and n = • 

III. Analytical calculation of stability zone of ZC1- 

DPLL WITH DESIGN LIMITATIONS 

The stability criterion of the non -ideal ZC 1 -DPLL can be 
obtained using the method outlined in the literature [7,8]. 
Define X(k) as a state vector given by 

X(k) = (a(k),P(k)), where a(k) and f)(k) are defined as 

a(k) = <p(k) and f3(k) = <p(k + 1) . If X* is a value of X such 
that X = G(X*). Then X* is called fixed point of G, and under 
certain conditions, the sequence [X(k) } will converge to the 

solution T, i.e., J™ 1 X(k) = X' . The Jacobian of the matrix is 
given by 

'0 1 1 

(12) 

J] P) 

Where, p = (1 - an) - aG^Kg cos [/?(&) + nfy\ and 

rj = an- G£K {1 + (n - l)a 2 } cos[a(£) + nfyr\ . 

Now, for frequency step input £ * 1 . At the steady state 

<p(fc + 2) = <p(fc + 1) = <p(k) = (p ss (say), gives the steady state 
phase error as 

. _! 2^(l + an)(g-l) 

= sin ; K^r 

^^{1 + a + (n-X)a 2 } 

For conventional DPLL, taking phase step input (£= 1) the 
steady state phase error ( q) ss ) is (zero) whereas for non- 
ideal DPLL (p ss = -nfy 

Now considering this for x* = q> ss , equation (12) reduces to 
(0 1 

G (X*) = 

[n s p., 

Where we define the following quantities, 

tl„ =an-^ (pr G, {1 + a + (n - l)a 2 ]f - (2^(1 + an){$ - \)f , 
p s = (1 - an) - a^GJl + a + (n - l)a 2 }f - (2tt(1 + an)(| - 1)) 2 
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7 = G'{X (k)} 



(13) 



Determining the eigen -values of the matrix (13) and applying 
the convergence condition [7], one gets the following stabil- 
ity condition 

2,0+^-1) < 1+an) j^Z±NE p! (14) 

G,{l + a + («-l)a 2 } G,{l + a + (n-l)a 2 ) V ' 

Equation (14) reduces to that of an ideal DPLL for a=l, r=0 
(i.e. G =1) and n=0 [7,8]. Thus keeping normalized input 
frequency (£) fixed, the phase locked condition could not be 



attained if the loop gain (K ) be more than 



V(l + a ) 2 +(2^(|-l)) 2 
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Again if keeping normalized input frequency (Q fixed, the 
phase locked condition could not be attained if the loop gain 

(K ) be less than 2x{i; - l)/G£ ■ 

IVA. Numerical Simulation Results 

The behaviors of the non-ideal ZC1- DPLL have been 
studied using the system equation (11). Fig. (3a) shows the 
bifurcation diagram of ideal ZC1-DPLL taking 1. Fig.s 
(3b), (3c) and (3d) show the bifurcation diagrams of non- 
ideal ZC1-DPLL for different values of design parameters 
(i.e. a and r) taking £ = 1 . From Fig.3b we have seen that for 
phase step input = 1) taking normalised delay parameter 
(a) =0.8 and normalised Sampling pulse width (r)=0 the phase 
locked condition could not be attained if the loop gain K Q 
more than 1 .72 and Fig. (3d) shows that for phase step input 
(E, = 1) taking a=0.8 and r=0.2 the phase locked condition 
could not be attained if the loop gain K Q more than 1 .92 and 
Fig.(3c) shows that for phase step input = 1) taking a=l 
and ?^0.2 the phase locked condition could not be attained if 
the loop gain (K Q ) greater than 2.2. The above results indicate 
that the upper limit of lock range increases due to increase of 
Sampling pulse width and decreases due to increase of loop 
time delay. Hence we may conclude that in a practical ZC1- 
DPLL circuit the upper limit of the lock range may remain 
almost unchanged as ideal ZC1-DPLL circuit for a suitable 
value of designing parameters a and r. It means that, although 
the presence of an additional delay reduces the stability zone 
of ZC1-DPLL, the finite pulse width increases the stability 
zone. Thus, a designer may make a trade-off between these 
two parameters to achieve an optimum loop operation. Fig. (5) 
shows this situation in the K Q - 1, space. 

Fig. (6) shows the fastest convergence of ZC 1 -DPLL with 
design limitations taking different values of designing 
parameters (i.e. a and r). It can be found that for an ideal ZC1 - 
DPLL, K=l should be taken for fastest convergence. For a 
ZC1-DPLL with finite width of sampling pulses and for a 
delayed ZC1-DPLL the fastest convergence to a steady state 
can be achieved for K Q >1 and K Q <1, respectively and for 
non-idealZCl-DPLL (i.e. considering both design limitations) 
taking a=0.9 and r=0.2 the fastest convergence to a steady 
state can be achieved for K =l and which is found in case of 
ideal ZC 1 -DPLL. It is also clear from the Fig.(6) that the effect 
on convergence time due to sampling pulse width and time 
delay oppose each other and for suitable values of parameters, 
a and r the overall effect may nullifies each other. 
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Fig.3 
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a Bifurcation diagram, plotting steady state phase error (p 
with different loop gain (K ) for an ideal ZC1-DPLL (i.e., 
r=0, n=l and «=0) at t, = 1 
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Fig. 3b Bifurcation diagram of (p ss (rad) with different loop gain (A" ) 
for an non-deal ZC1-DPLL having inherent time-delay o=0.8 
(with r=0 and «=1 ) at % = 1 
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Fig. 3c Bifurcation diagram of (p ss (rad) with different loop gain (S" {) ) 
for an non-deal ZC1-DPLL having finite pulse width (r=0.2); 
(with a=l and n=l) at § = 1 
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Fig. 3d Bifurcation diagram of (p ss (rad) with different loop gain (K lt ) 
for an non-deal ZC1-DPLL (i.e. n=l) having inherent time-delay 
(a=0.8) and finite pulse width (r=0.2) at I = 1 



IVB. Time Series Analysis 

The occurrence of chaos in a dynamical system can be 
qualitatively measured by computing average Lyapunov 
exponent. Lyapunov exponents quantify the exponential 
divergence of initially close state-space trajectories and 
estimate the amount of chaos in a system. A positive 
Lyapunov exponent indicates chaos and in this situation the 
system will be very sensitive to initial conditions. The 
techniques of calculation of Lyapunov exponent from time 
series data of a system variable is well documented in 
literature [9,10]. We calculate the Lyapunov exponent of ideal 
ZC 1 -DPLL and non -ideal ZC 1 -DPLL from time series data of 
phase errors with different loop gain values K Q (for phase 
step input, i.e. Jj= 1). Fig. (4a), (4b), (4c) and (4d) show the 
Lyapunov Exponent of ZC1-DPLL for different values of 
designing parameters (i.e. a and r) with £ = 1. Fig. (4b) shows 
that taking a=0.8 and r=0 the average Lyapunov exponent is 
negative as the controlled parameter K Q is varied from to 
1.72 and Fig. (4d) shows that for a=0.8 and r=0.2 the average 
Lyapunov exponent is negative as the controlled parameter 
K Q is varied from to 1 .92 and the same result also has been 
found analytically. Fig. (4c) shows that taking loop delay 
parameter a=\ and normalized sampling pulse width r=0.2 
the average Lyapunov exponent is negative as the controlled 
parameter K Q is varied from to 2.2. From Fig. (4b) and Fig. 
(4d) we conclude that the reduction of steady state phase 
error zone due to loop time delay is more or less compensated 
due to increase of sampling pulse width. Here we may 
conclude that the results found related to steady state zone 
using Lyapunov exponent agree with the results found from 
bifurcation diagrams and analytical predictions. 
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Fig. 4d Plotting of average Lyapunov Exponent (A, ) with 
different values of loop gain K for (a) ideal ZC1-DPLL (b) a=0.8, 
r=0 (c) 0=1, r=0,2 (d) 0=0.8, r=0.2 = 1) 
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Fig. 5 Stability zone for ideal ZC1-DPLL and non-ideal DPLL with 
values of design limitations parameter a=0.8, r=0.2 
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Fig. 6 Plotting of convergence time (in number of cycle) with 
different values of loop gain A" for ideal ZC1-DPLL and ZC1-DPLL 
with finite width of sampling pulses 



V. Conclusions 

This paper reports the effect of loop time delay and finite 
sampling pulse width on the dynamical behavior of a ZC1- 
DPLL. For an ideal ZC 1 -DPLL one considers that there would 
be no additional loop time delay and also the width of 
sampling pulses are narrow enough. But in a practical system 
there is always an additional time delay in the loop and also 
the sampling pulses are not instantaneous but they have a 
finite width. These two effects should have to be considered 
to explore the complete dynamical behavior of a ZC1-DPLL. 
The stability of the loop has been found to be affected by the 
loop time delay and finite width of sampling pulse. Further 
using nonlinear dynamical loops like bifurcation diagram 
(local) and time series analysis it has been shown that the 
loop dynamics shows period doubling bifurcation and chaos 
for the variation of loop time delay and sampling pulse width. 
The authors believe that the present study would help the 
designers to design optimum practical DPLL where the 
practical design limitations (such as loop time delay and Af ) 
have been considered. 
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